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Using post-Newtonian equations of motion for fluid bodies valid to the second post-Newtonian 
order, we derive the equations of motion for binary systems with finite-sized, non-spinning but 
arbitrarily shaped bodies. In particular we study the contributions of the internal structure of the 
bodies (such as self-gravity) that would diverge if the size of the bodies were to shrink to zero. Using 
a set of virial relations accurate to the first post-Newtonian order that reflect the stationarity of 
each body, and redefining the masses to include IPN and 2PN self-gravity terms, we demonstrate 



Othe complete cancellation of a class of potentially divergent, structure-dependent terms that scale as 
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and s~^^'^ , where s is the characteristic size of the bodies. This is further evidence of the Strong 
' Equivalence Principle, and supports the use of post-Newtonian approximations to derive equations 

of motion for strong-field bodies such as neutron stars and black holes. This extends earlier work 
done by Kopeikin. 



I. INTRODUCTION AND SUMMARY 

The principle of equivalence is the cornerstone of Einstein's general theory of relativity. Part of this principle, 
^ known as the weak equivalence principle (WEP) states that test bodies fall in a gravitational field with the same 
I acceleration, irrespective of their internal structure or composition. By test body, one means a body whose internal 

structure is governed by non-gravitational interactions, and whose size is small compared to inhomogeneities in external 
■ gravitational fields (suppression of tidal couplings). WEP, together with postulates of local Lorentz invariance and 
^ local position invariance of non-gravitational physics in freely falling frames, comprises the Einstein Equivalence 
£2 Principle (EEP), which is the foundation of metric gravity, or of the idea that gravity is really geometry. 

] There is a stronger version of WEP, which proposes that all "test" bodies should fall with the same acceleration, 
. even bodies that are self-gravitating, such as stars, planets, or black holes. Here, by test body, one means a body that 
is large and massive enough to have measurable self-gravity, yet is small enough that tidal interactions can be ignored 
(one generally ignores spin interactions as well). This version of WEP is a part of the Strong Equivalence Principle 
(SEP), which includes a postulate of a kind of local Lorentz invariance and position invariance of gravity itself. While 
\ every metric theory of gravity obeys EEP, almost no metric theory obeys SEP. Scalar-tensor theories of gravity, such 
as the Brans-Dicke theory and its extensions, violate SEP. Indeed, SEP is generically violated in theories of gravity 
that introduce long-range fields, in addition to the spacetime metric, that mediate how the metric is generated by 
matter. In order to preserve EEP, such fields do not couple directly to matter. Because general relativity contains 
^ ' one and only one long-range gravitational field - the metric itself - it has no mechanism for violating SEP. Far from 
any gravitating system, the metric can always be made suitably, if approximately, Minkowskian, and so there is no 
obvious mechanism, other than tidal interactions, for the external universe, or for any nearby "spectator" body to 
influence the internal structure or dynamics of the system (for further details of this argument, see [l[ and Sec. 3.3 
of Q)- Thus gravity in GR is independent of the velocity of the system relative to some external frame, the effective 
constant of gravity, G, is a true constant, and self-gravitating non-spinning bodies move as if they were test particles 
(for further discussion of SEP, see Section 3.3 of (^f; for an alternative discussion, see Q). 

But these are general, qualitative arguments. This paper addresses the question: how explicitly does general 
relativity manage to satisfy SEP for self-gravitating bodies, specificially to the second order in a post-Newtonian 
expansion? 

At the first post-Newtonian (IPN) order of approximation, that is, at order e ~ (w/c)^ ^ Gm/rc? beyond Newtonian 
gravity, GR has been shown to obey SEP explicitly. This is seen most graphically within the parametrized post- 
Newtonian (PPN) framework, which characterizes the post-Newtonian limit of a wide range of metric theories of 
gravity using 10 arbitrary parameters (see Chapter 4 of [4I for a review). Some of these parameters or combinations of 
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them measure whether the theory has "preferred-frame" effects at post-Newtonian order (violations of local Lorentz 
invariance) , some measure whether the locally measured gravitational constant can depend on the presence of nearby 
matter (violations of local position invariance), and some measure whether self-gravitating bodies violate WEP, a 
phenomenon known as the Nordtvedt effect. In GR all such offending parameters or combinations of parameters 
vanish. Furthermore, there is experimental evidence that they vanish from a variety of tests of post-Newtonian 
gravity, spanning lunar laser ranging, geophysical measurements and binary pulsar data (see Q for a review). 

But what about beyond post-Newtonian order? Does the motion of a pair of finite sized, gravitationally bound, 
non-spinning bodies depend on their internal structure, apart from tidal interactions? There is existing theoretical 
evidence that they do not. Breuer and Rudolf Q showed that the relative "Newtonian-like" acceleration of two well 
separated bodies, momentarily at rest, was independent of their internal structure, to second post-Newtonian (2PN) 
order, or order beyond Newtonian gravity; put differently, they proved the absence of the Nordtvedt effect to that 
order. Kopeikin [6] derived the equations of motion to 2PN order for self-gravitating bodies with an equation of state 
p{p), under the assumptions that they were spherically symmetric in a suitable comoving frame, and had no internal 
fluid motions, and also found that the SEP held (see also [3])- 

This question is of more than academic interest. The leading candidate sources for gravitational radiation that may 
be detected in the near future by laser interferometers are binary systems containing neutron stars or black holes. 
The inspiral part of their evolution can only be described accurately by the post-Newtonian approximation, which, 
by its very nature, assumes that gravitational fields are weak everywhere in spacetime. No attempt has ever been 
made to derive the motion to high orders in a PN expansion with either strong-gravity neutron stars or black holes. 
The closest one has come is to treat the sources as distributions in spacetime (delta functions) and to employ suitable 
regularization techniques to control the infinities d, Q , or to use a variant of the surface integral approach pioneered 
by Einstein, Infeld and Hoffmann [lo| . Therefore an explicit demonstration that the internal structure is irrelevant 
for the binary motion, even if only at 2PN order, would be valuable. 

In this paper we derive the equations of motion of two arbitrarily shaped, non-spinning, gravitationally bound 
bodies, through 2PN order. Beginning with the 2PN hydrodynamic equations of motion derived in Paper I of this 
series 111 , we calculate the acceleration of the center of "baryonic" mass of a given body in the field of a companion 
body. Expanding all variables about the centers of mass of each body we track all terms that scale with the sizes s of 
the bodies as s~^, and s^^/^; such terms represent contributions due to self-gravitational binding energy, and a kind 



of gravitomagnetic internal energy, respectively. (In Paper II |12| we ignored all finite-size effects in deriving binary 
equations of motion.) We then employ virial relations that refiect the equilibrium state of each body, leading to a 
cancellation of many terms. Finally we renormalize the masses to include both rest mass and gravitational binding 
energy, by defining the total mass for body A to be 

Ma = niA + ^CIa + OiniAe^) , (1.1) 
where uia = J a P*dPx is the conserved baryonic mass of body A and 

nA=-l [ [ J^d'M'x' (1.2) 
2JaJa\^-^'\ 

is the gravitational self-energy as measured in a local, freely falling frame x", momentarily comoving with body A. 
The result is that all self-gravity terms that scale as s~^ or s~^/^ cancel, leaving equations of motion at 2PN order 
for body 1 that depend on the rcnormalized masses Ma, given by 



r 

where r = |xi — X2I, 71* = {xi — X2f jr, and where 
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-4(vi • n)(vi • V2) + 4(v2 • n)(vi • V2) + -(v2 • n)^ 



(vi • n) + (V2 • n) 



(1.4b) 



The equation of motion of body 2 can be found by interchanging 1^2. These agree completely with other results 
for the 2PN equations of motion [1 S d, [H [ll [11, HI . 

The remainder of this paper provides details. In Sec. |lT] we lay out the basic equations and assumptions that 
underlie our method. Section [IIII focuses on the first post-Newtonian approximation and verifies SEP to that order, 
while Sec. llVl extends this to second post-Newtonian order. Concluding remarks are made in Sec. |Vl Selected detailed 
calculations are relegated to a series of Appendices. 

II. BASIC EQUATIONS AND ASSUMPTIONS 

We consider the motion of a binary system made up of two fluid balls of characteristic mass m and size s. We do 
not assume that the bodies are spherically symmetric. They could be distorted, either because of rotation or because 
of tidal interactions (although we will not take either spin or tidal effects explictly into account). The only symmetry 
we will impose on the bodies is that each be symmetric on reflection through its center of mass. In practice this 
means that we will ignore any odd-£ angular harmonics in its density distribution; equivalently we will assume that 
the integral of any odd number of purely internal vectors (such as position or velocity relative to the center of mass) 
over the body vanishes. 

We will treat the material making up the bodies as a perfect fluid, with the pressure required for equilibrium 
provided by random internal fluid motions. Thus we will break up the velocity of each fluid element in the A-th 
body according to v = v^i -I- v, where v^i will be a suitably defined center-of-mass velocity, and v could contain both 
random thermal- like velocities as well as bulk internal or rotational velocities. We set the formal pressure equal to 
zero (j) = 0). We also characterize the fiuid by the so-called "conserved" baryon mass density, p*, given by 



,0 



mn\/—gu 



(2.1) 



where m is the rest mass per baryon, n is the baryon number density, is the four-velocity, and g = det((7^,y) is 
the determinant of the spacetime metric g^^,. Henceforth, we use units in which G = c = 1; Greek indices range 
over spacetime values 0, 1, 2, 3, while Latin indices range over spatial values 1, 2, 3. From the conservation of baryon 
number (expressed covariantly as Vainu") = 0), p* satisfies the exact continuity equation 



dp* 

^+V.(p*v).0, 



(2.2) 



where v" = u°' /u^ — (l,v) is the fluid's coordinate velocity, and spatial gradients and dot products use a Cartesian 
metric. In terms of p* , the stress-energy tensor takes the form 



T 



Q/3 



p-{~9) 



(2.3) 



We define the baryon rest mass, center of baryonic mass, velocity and acceleration of each body by the formulae 



ruA 



VA 



SiA = 



p*d^x , 



(l/m^) / p*^d^x, 

I A 



dxA/dt — (I/ttia) / p*vd x, 

J A 

dw a/ dt — {l/niA) / p*aLd'^x, 

J A 



(2.4) 
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where we have used the general fact, imphed by the equation of continuity for p*, that 

^ J p*{ty)f{t,^,^')d^x' ^ J p*(i,x')(^^+v'.V'^/(t,x,x')dV. (2.5) 

The definitions of mass and center of mass are not unique, of course (for a review, see Q)- For example, we could 
have used an effective density that included internal kinetic and gravitational potential energies, as in Sec. 6.2 of 
0, for example. As it turns out, our assumption that only even £ multipoles of the internal density distribution 
matter guarantees that all such alternative definitions actually coincide with our baryonic definition, by symmetry. 
Nevertheless, one could generalize our assumptions and consider the effect of different center-of-mass definitions, but 
that would be beyond the scope of this paper. 

To 2PN order, the equations of hydrodynamics have the form [II, Eq. (2.23), (2.24)] 



dv'/dt = U'' + a'pj^ + a'2PN , (2.6) 



where 



J-PN 



'■2PN 



1 (3) 

+w^t;'=(2$f - 4$f'^ + 2Pf - AP^'^^) + i;'(3$2 " " ^ ^ +W''U''') 

-Aq>Y - 4P;^) + ^S'Xw*) + ^J:^'{v^U) ~ ^(y^V^) - ^S'^^i) - 6[/$i' - 2$iC/'^ 
8 2 2 

_4$5^'[/J + 8t/^T/^'' + 4F'f7 + 2t{v^v^) + 4C/$2 + 4$2C^'' + SC/^L/'* - S^*($2) 

3 .. .. . 1 .. 3 .. . 1 .. . (3) 1 (4) 

+ -I]''(L/2) - 2C/X'* - 2X[/'' - WV - -E'X^) + -X-. - -Xk' + 2 X* +— Y^' 
2 2 4 2 24 

+AVi - 8$^ - QG'l - 4G^* + 8G3 + 8G4 - 4Gg + 16GJ, - APl^U^^ - , (2.7b) 

where commas denote partial derivatives, overdots and the notation (n) over functions denote partial time derivatives 
d/dt, parentheses and square brackets surrounding indices denote symmetrization and antisymmetrization, respec- 
tively, and repeated spatial indices are summed. 

Definitions of the potentials appear in Appendix [XJ all are defined using p*. For example, U is the Newtonian 
potential, defined by 

where is a constant-time hypersurface extending to the boundary of the near-zone of the system. 
Our task is then to calculate the acceleration of, say, body #1 through 2PN order, using 

a\ = (1/mi) / p*{dv'/dt)d^x. (2.9) 

At Newtonian order, the result is simple, namely 

a{ ^ -{l/m,) jj^p*p*'^^^d^xd^x' 

+{i/mi) J^p*d^x J^p*'d^x' jvii + ix- x'yv'/^ + ^(^ - x y'^^Y''^ + . . . | , (2.10) 

where we split the position and velocity of a given point inside each body A according to 

= v\+v\ (2.11) 
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and expand the potential due to body #2 in a Taylor series about the centers of mass of the two bodies. The use of 
multi-indices denotes products of vectorial objects; for example = x'^x^ , V*-' = V'V^ , x^ = . . . a;*™ , and so on. 

The first term in Eq. (|2.10p vanishes by symmetry (Newton's third law). The first term in the second line 
corresponds to the normal "point-mass" Newtonian acceleration, —m2n^/r^; the second term vanishes by the definition 
of baryonic center of mass, and the remaining terms are standard Newtonian multipole coupling terms. They depend 
on the size of the bodies as s" with n > 2. In the limit that the size of the bodies is negligible relative to their 
separation, the latter terms are vanishingly small; these are the kinds of terms that we will ignore throughout. 



III. EQUATIONS OF MOTION AT IPN ORDER 



At higher PN order, there is now the possibility of correction terms proportional to inverse powers of s. To illustrate 
this, we evaluate two terms from Eq. (j2.7ap explicitly: 



— / p*v^U^'d^x 
mi J I 



— [ p*(vl + 2vi-v + v^)d^x 
mi Ji 

p*'{x-x')\^ , m2x\2 m2x\6'^ -3n'n^) 
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2Ti m2x\ 
mi 



(3.1) 



and 
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mi Ji 



p*UU''d^x = 
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il^^ m2Xi2 _|_ 2 ^-'^ "^2''2^i2 m'^x\2 



mi r" 



mi 



0{s), 



(3.2) 



where 71, fl^ , fli and T-C^ are defined in Appendix [Bl Note that, because ^ m/s for a body in equilibrium, the 
"kinetic energy" 7^ scales as 1/s, as do the "gravitational potential energy" quantities ili and fl^ ■ The "gravitomag- 
netic" style quantity Tii scales as v/s"^ ~ l/s^/^. We have used our parity assumption to eliminate all terms that 
involve an integral over an odd number of internal vectors on a given body. Also, we have not kept terms that scale 
as positive powers of s. One example is a term from Eq. (|3.ip proportional to xv ^ s^/^; for rotating bodies it gives 
spin-dependent effects, which we treated in Papers III and IV 17, 18]. Here we ignore spin. 

In the combination of IPN terms 4V^ + \X'^ in Eq. (|2.7ap . the time derivatives generate accelerations inside the 
potentials. To the order needed for a IPN calculation, we must therefore substitute the Newtonian hydrodynamic 
equations for those accelerations and carry out the same procedures for the integrals as described above. When we 
go to 2PN order, we will have to insert the IPN hydrodynamic equations. 

Carrying out these procedures for all the terms in Eq. p.7ap . and dropping terms scaling as positive powers of s, 
we obtain for the equation of motion of body 1 , 
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(3.3) 



Note that, apart from the "point" mass terms, only terms scaling as s^^ and s^^/^ arise. 

We now assume that each body is in equilibrium, implying that it is stationary, or possibly periodic on an internal 
dynamical timescale. As a consequence, any time derivative of the moment of inertia tensor, = J p*x'^x^d^x, can 
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be set to zero or can be averaged to zero. This results in a set of virial relations derived in Appendix |B1 which imply 
for each body that 2T^ + fi^ = 0, 2Ta + = 0, and AHa'^ — = 0. These eliminate most terms dependent 
upon the structure of the bodies, leaving only the term — nV~^(472 + ^^^2)- This cannot be eliminated by a virial 
relation; however we can add an arbitrary amount of 272 + 5^2 = to it to put it into the form 



(4 - 2a)r2 + (- - a)fl2 



(3.4) 



Despite these terms, we can make the IPN equations independent of any 0{s ^) self terms by redefining the masses 
of each of the bodies to be. 



= + (4 - 2a)TA + 



a n. 



(3.5) 



resulting in Eq. (|1.4ap . Note that the redefinition of the masses from tua to Ma in the IPN terms will affect the 
equations of motion only at 2PN order. Thus we have verified the SEP to IPN order. 

The choice a = 3/2 gives a redefined mass Ma = ruA -\-Ta + ^^^7 which naturally represents the total (baryonic plus 
kinetic plus gravitational) mass of the body, while the choice a — 2 gives Ma — itia + ^^a, which is the same thing, 
after applying the virial relation Ta — —51^/2, and represents rest mass plus gravitational binding energy. Either 
definition, when applied to an isolated gravitating system, gives the total system mass, as measured using Keplerian 
orbits far from the system. At IPN order, the choice of a is completely arbitrary, but we will see that, to satisfy SEP 
at 2PN order, we must choose a = 2. 



IV. EQUATIONS OF MOTION AT 2PN ORDER 



We now apply the same methods at 2PN order. We restrict attention to terms analogous to those that arose at 
IPN order, namely terms that scale as s~^, s~^/^, as well as "point-mass" terms. The most straightforward terms 
to evaluate are those that involve two-body potentials and superpotentials, such as ?7, V, $1, A, and their various 
generalizations, such as V2, X'\ Y, etc (see Appendix [K\ for definitions). We use the splitting of position and velocity 
as in Eq. (|2.1ip . combined with the scaling v ^ s~^^^ and x ~ s. It is a simple matter to expand and sort the terms 
using algebraic software. To illustrate the results, we cite a few simple examples from Eq. (j2.7bp : 

— / p*v'v'v''V''''d'^x = -A^vlivi ■ n)(vi • V2) - — ^ f wjw^n'^'r/'^ + (vi • n)viT^' + (vi • V2)n'T^- 
mi Jl mi V 

+ — ivlvllCi + v{vlre') , (4.1a) 

TOl 



ii Jl mir^ 



(4.1b) 



Throughout, we apply the Newtonian virial relations, 2T^ + fi^ = and 47i^"''' — 3/C^ = to all 2PN terms that 
scale as and s~^/^. In some terms, an additional virial relation is required, involving d^V^/df^. An example is 
the term in Eq. (|2.7bp involving d'^Y'^/dt'^, where Y is the superduperpotential (see Appendix [XI for definitions). 
Splitting l"(x) into two terms, one from body 1 and one from body 2, one can show, using a Taylor expansion of the 
second term about the center of mass of body 2, that, inside body 1, y is given by 

y(x)= / p*V-xf dV+TO2|x-X2p + i/fV^'=|x-X2|3 + -^/f'"V^'='™|x-X2|3 + ... . (4.2) 

J I 2 12 

Partial time derivatives of the terms involving the moments of body 2 either will give zero, because of our virial 
assumption that the moments are effectively constant in time, or will give terms proportional to the moments them- 
selves, which therefore scale as positive powers of s. Thus in this case, only the contribution from body 1 and the 
point mass term from body 2 will contribute; no self-terms from body 2 will arise. 

Next in complexity are the so-called triangle potentials, such as P2 = P{U''^U'^), and the potentials Ga {a = 1..6) 
and G\, which depend on the field point and on two source points, hence the name triangle potentials. Luckily these 
can all be written in analytic form using the "triangle" function Q{xAB), which satisfies the differential equation 
'V^Q{xAB) = — (|x — x^||x — xb|)^^, and which is given explicitly by 

g{xAB) EE ~lnA{xAB) + l, 
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A{xAB) = |x-xa| + |x-xb| + |xa-xb|. (4.3) 
For example, the potential G\ = P{U'^) can be expressed in the form 

Gi^J2 f I P*AP*Bd''xAd?XBv\v'gV\V'Bg{xAB) . (4.4) 



A,B 



Consequently, terms involving triangle potentials can be evaluated with ease (using algebraic software, to be sure) 
using the same splitting and expansion procedure as before. For example. 
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m 



TtT-i 

+ [3nXvi • V2) - 12n^(n • vi)(n • V2) + 6vl{n ■ V2) + 3w^(n • vi)] — 

r 



^(n • V2){vin, ~ vi^i^) ~ 6^{vln2 - viv^ni") . (4.5) 



The most difficult term to evaluate is the final term in Eq. (j2.7bp . — 4iJ'\ This involves the so-called "quadrangle" 
potential, since H — P{U'^^ P2) is a function of the field point and three source points. Unfortunately there is no 
known analytic formula similar to the function Q that can be employed to simplify this potential or to expand it about 
the bodies' centers of mass using algebraic software. Instead an alternative and cumbersome method must be used; 
this is detailed in Appendix iDl The result for this term is 

_ A / p*iJ.M3x = _!!^(8mim2 + m2) 

^{ (5mi — TO2)n* — - + (5mi + 1712)11^ — 2(10mi — m2)n^n^n''^— 

mi mi mi 

+2n'n2 + 2{mi+3m2)n^— ~2{mi+6m2)n''n'n^— )■ . (4.6) 

TO2 7712 J 

In the limit of spherically symmetric bodies, where iT^ = HaS^-' /i, all the contributions cancel, and the result 
agrees with that of 

We also must insert the IPN hydrodynamical equations of motion into the accelerations that appear in the IPN 
terms AV'' + \X''^ in Eq. (|2.7ap . and evaluate those 2PN terms. They involve only 2-body potentials, and thus are 
straightforward to evaluate. 

Combining all the terms that scale as and s'~^/'^, and displaying only those terms, we obtain the result 
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-71^1 (4vi • V2 + 20^) ^ + (37;^ + 15 
+ («i - V2r\ f 2n . vi - . V2') ^ + Av{n^^\ 

IV 2 / 7712 7711 J 

+ ( 4(n • Vi)u^ - (n • Vi)7;;? - (n • V2)7; 1 ) — 

V / TTtl 



- -vM^-^ -v{v\v[ { -^f ' - 3/:f' - erf-^' 

TTii 2 m2 mi V 2 



(4.7) 



where ri^*^', which scales as s ^ , and fC^^^ and C^^^^ , which scale as s are defined in Appendix [Bl 

We now must return to the IPN equations of motion including and s~^/^ terms, Eq. (|3.3p . and employ virial 
relations that are correct to IPN order. As we saw in Sec lIIIi all self terms vanish at IPN order, except for the term 



8 



proportional to 472 + 5/2 fi2- The residual term there could be absorbed into a redefinition of the masses. We must 
now consider the application of virial relations and the mass redefinition at higher PN order. 
We first replace the term —{n'^/r^){4:T2 + 5/2 il2) in Eq. (|3.3p by the equivalent term 



(4 - 2a)r2 + (- - a)n2 



(4.8) 



and apply the IPN corrected virial relations of Eqs. (jB5p only to the second piece of this expression, as well as to all 
the other and s~^^^ terms in Eq. (|3.3p . The other piece of (|4.8p will be absorbed into a redefinition of the mass of 
body 2. At the same time, we redefine all the masses in the point-mass IPN terms using Eqs. ()3.5|) : for those terms, 
to the order of approximation needed, the Newtonian virial relations may be used to simplify the renormalization to 
Ma — TUA + 1/2 ^Ia- The s^^ and terms generated by these substitutions cancel all the s~^/^ terms and almost 

all the s^^ terms in Eq. (|4.7p . The uncancelled terms at IPN and 2PN order, combined with the Newtonian 
acceleration, give 



ai(self) 



TO2 + (4 - 2a)r2 



-{a - 2)vln2 - i(21 - 10a)^r!2 " ^(13 - 6a)viv^ni'' 



(4.9) 



Note that the choice a — 2 leaves the coefficient of the "Newtonian" acceleration as 

1 



m2 



-^2 

2 



-^n2 

2 r 



V2'^2"2 



= m2 - 

= M2.. 



(4.10) 



where 



^]n2 

r 



^2^2 



(4.11) 



is the gravitational binding energy as calculated in the local, comoving inertial frame of body 2 (see Appendix [C] for 
derivation). The quantity M2 is precisely the total mass, comprising the baryonic mass m2 plus the locally measured 
gravitational binding energy ^Ct2. The additional terms in (j4.1ip are simply the corrections arising from the 
transformation from the local, comoving inertial frame to our global coordinate frame. 

Thus, when all masses are written in terms of the new Ma, no terms survive in the equations of motion, leaving 
only the "point"-mass terms, given by Eq. (|1.4bp . 



V. CONCLUSIONS 



We have studied the motion of finite sized, self-gravitationally bound, arbitarily shaped bodies at second post- 
Newtonian order, and demonstrated that self-energy effects that scale with size of the body as and s"^/^ cancel 
precisely when a suitable renormalization of the masses is carried out. Since the bodies are finite, that renormalization 
is finite (in contrast to other, singular renormalization techniques), and consists of redefining the masses, to the PN 
order needed, as the sum of baryonic rest mass and gravitational binding energy, as measured in the locally comoving 
inertial frame of each body. We emphasize that we have made no effort to define masses or centers of mass in a 
covariant manner; all calculations are carried out in the harmonic coordinates that are built into our approach. We 
are currently investigating 2PN terms with other scalings, such as s~^/^, and so on, to see if they also cancel. 
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APPENDIX A: POTENTIALS APPEARING IN THE 2PN EQUATIONS OF MOTION 

The potentials that appear in the equations of motion are all Poisson-like potentials and their generalizations, 
namely a superpotential X and a superduperpotential Y , integrated over a constant time hypersurface AA that 
extends to the boundary of the near zone of the system. In the case of integrands of non-compact support, we discard 
all terms that depend on the radius of the near-zone, 7^; see Paper I [ll[ for justification. All potentials are defined 
in terms of the conserved baryon mass density p*: 

JM |X - X I 

X{f) ^ f p*(t,x')/(t,x')|x-x'|dV, 

JM 

Yif) EE f p*(t,x')/(i,x')|x-x'|W. (Al) 

JM 

The specific potentials used in the IPN and 2PN equations of motion are then given by 

$1 = $2 = S(C/), ^ = ^(1), 

= ^vV), $^ = £(1/*), Y = Yil), 
X' = Xi = Xiv^), X2 = X{U), 

Gi = P(C>2), G2 = P[Uty), 

G3 EE ~P{U^^V^), Gi = P{V''W'''), 

G5 ^ -P{V''U'''), Ge = P(C/'''^$i^), 

G^ = P{U'''V''^') + ^PiU-'U), H = P{U^'^Pi^). (A2) 
We refer the reader to Ref. Appendix C, for further discussion of the triangle and quadrangle potentials. 



APPENDIX B: VIRIAL THEOREMS 

1. Newtonian virial relations 

We assume that our bodies are in equilibrium, so that they are either stationary, or at worst periodic on an internal 
dynamical timescale. This implies, among other things, that any time derivative of the moment of inertia tensor, 
= J p*x^x^ (Px, is either zero or averages to zero. This will give us several virial relations that will simplify our 
equations. Considering body 1, for example, we have that 

i/f = [ p*v<^'x^U^x, (Bla) 

2 Ji 

ij;^' = [ p*{v'v3 +x^'a3^)d^x, (Bib) 

2 Ji 

I'ii' ^ f p*i3v'^'a^^ +x^'d^^)d^x, (Blc) 

2 Ji 

where, by virtue of the fact that p*x'^d^x = p*v^d^x = 0, we can drop the bars on the accelerations. Substituting 
the Newtonian equation of motion for a* gives 

= ^ p*v'v'd'x - JJ^ p*p*'^-^^^^d^xd^x' + ^ p*x^'U^U''x , 
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= 2T;' + + 0{s^) , (B2a) 



where 



1 Ji 



2-/^ = -/ /pvm3 + -3::-^ — ^n^^zT^^^ — ^H^-'^v 



4Hf ^ - 3IC\' + 0(si/2) ^ (B2b) 



1 Ji 

'=? s / / pV" " ''}""' 'if ~ d'xd'x' . (B3) 



1 Jl 



and where we expanded the potential U2 in the same manner as in Eq. p.lOp . yielding only terms that scale as 
positive powers of s. (In Paper III [l7j . the s^/^ terms in Eq. (|B2bp produced some spin-orbit terms that contributed 
to the IPN spin-orbit equations of motion.) Recall that, in these definitions, {x — x'Y^--- = {x — x')''{x — x'Y .... 
Setting the derivatives of to zero, we obtain the and s~^/^ virial relations for body 1, 

2T^'+n\' - 0, 

mf'^~3IC\' = 0. (B4) 



2. Post-Newtonian virial relations 



Working at 2PN order in the equations of motion requires us to obtain our virial relations correct to IPN order. 
Accordingly we must now substitute the IPN hydrodynamical equations into Eqs. (jBip . The resulting IPN corrections 
will have terms with a variety of scalings, from purely internal terms scaling as s^^ in and s^^/^ in to terms 
arising from expansion of the potentials due to body 2, also with a variety of scalings. Here we focus only on 
contributions scaling as s^^ and s~^/^. A straightforward calculation then yields the IPN virial relations, 

- 2T;' + - vlnY ~ '^-v\v\'^'l''' - h^^'l , (BSa) 



. + 4^(z;i - 2^2)'=n('f7f - g^^^^f^f n'^ + 3^.;^'^^' , (B5b) 

where v = vi — V2, and 



. //p>"'-^--^''""-'„'''-''''d3A'. 

= / \ p*p*'''^l^—^d^xd^x' . (B6) 



|x - x\ 

The virial relations for body 2 can be obtained from these by the interchange 1^2, with n 
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APPENDIX C: RENORMALIZED MASS IN BODY'S REST FRAME 

We have renormaUzcd the masses of the bodies by defining the total mass to be a sum of baryonic mass and 
gravitational binding energy, namely Ma = mA + modulo 2PN corrections. However, those 2PN corrections 
are of order me^ ~ m{m/s)'^, and thus scale as s~^. As we are only looking at terms, we have not kept those 
corrections. However, there is an additional, 2PN order, correction to the total mass that must be taken into 
account. 

The baryonic mass of each body is a scalar, frame-invariant quantity, but the gravitational binding energy is not, 
since it depends on the size of the body. The correct, frame-invariant definition of our total mass must therefore be 
that mass as measured in a local inertial frame momentarily comoving with the body. Thus we define 

Ma = mA + ^CIa + 6Ma , (CI) 

where 6Ma denotes the 2PN, corrections, which we are ignoring, and 



Iff p*p 



A*' 



where hats denote spatial variables defined on a constant time hypersurface in the comoving frame. Now, the quantity 
p*(Px is invariant, but |x — x'| is not. In Paper II, Appendix B, we showed that the transformation between spatial 
coordinates in the comoving frame of body A and our global harmonic coordinates takes the form 

x'^x\ + x^ {6] + e{A) + B)kx'')] + Oie") , (C3) 

where v4*- and S^j, are functions of the basis transformation from the global frame to the comoving inertial frame, 

e^^{Kl+B%x'')es., (C4) 

where — x'^ — x\. The coefficients A|^ correspond to boosts and coordinate rescalings, and the B'^^, make the frame 
freely falling. 

Substiting Eq. (|C3p into the definition of fii, we obtain 

f7i = f7i - tA'^ifl + 0{e^)Vli , (C5) 

where the term arising from i?*^. produces integrals over body 1 involving an odd number of vectors, which therefore 
vanish. The coefficients Aij arise from a boost to the velocity v\ of body 1, combined with a rescaling to an orthonormal 
basis in the field of body 2. To the first PN order needed, it is simple to show that 

A^ ^-^vlvi-'^S'' . (C6) 



The result is, to the required order. 



ni = hi + l-v\v'M + —Cii . (C7) 

2 r 



Therefore, to 2PN order, and with ,s ^ scaling, the total mass of each body is given by 

Ma - mA + ^^A + 0{s-^) 



= mA + ^flA - \v>'A^t - I'-^^A + Ois-') , (C8) 
where we drop the hats on the 2PN terms. 

APPENDIX D: THE QUADRANGLE POTENTIAL H 

We now turn to evaluation of the term involving the quadrangle potential H — P{U'^^ P2 )- This potential can be 
written in the form 
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= E / P^ViiV^^d^a;^ / p*By'Bd'xB [ ph^'cd^xcHixA; BC) , (Dl) 
at^^Ja Jb Jc 



ABC 



where the function Ti. of four field points is defined by 



Unfortunately, there appears to be no closed-form analytic expression for Ti similar to Eq. (|4.3|) for Q. Instead, we use 
the first definition of H in Eq. (|D1[) , and integrate H'^ over the mass of body 1 . After integrating once by parts, showing 
that the surface term at the boundary of M can be discarded, and using the fact that = ^^'^ - \UU'^ - 
we obtain 

[p*H''d^x^^f U{''U'''Pt£x + ^ f U{'U'''['P'2^ -UU^'' -2J:{U''')]£x, (D3) 
Ji 47r 8n 

where Ui is the Newtonian potential due to body 1 only. The first term in Eq. (jD3P can be expanded into the form 

Termi = i- |^ U^^^ (U^ + Utma,^ + ^^Si'i + ^22)]^'^ > (D4) 

where the subscripts denote the contributions from the various bodies. 

The second term in Eq. (jPSp can be simplified by integrating by parts, leading to 

Term2 = 7 / p*iU'<^2 + U^'^ - U'^U'')dPx + \ ( p*U{'^2d^x [ U(U^^T.{U'^)dPx , (D5) 

4 Ji ^ Jm 47r Jj^ 

where the first integral is only over body 1. 

To handle the integrals of noncompact support integrands, we split the domain M into three regions, a region A^i of 
radius TZi surrounding body 1 and a region M2 of radius TZ2 surrounding body 2, and the remainder, M. — Mi~ M2- 
To carry out the integrals, we will need suitable forms for the various potentials, J7, ^'2(11)' -^2(12)' ^"^^ ^''^^ 
appropriate regions. 

This will be aided by a general expansion of the function Q{ABC) in powers of tabItac, where points A and B 
are assumed to lie inside one body, and point C is inside the other body, so that tab ~ s <C tag- Straightforward 
methods lead to the expansion 



^ oc 

Q{ABC) = - InrAc + 1 - ln2 + - ^ 



2 ^ (m + 1)! 

m=0 ^ ' 



^UuABr^f f— ) + -^AuABr^'^T' f— ) j • (D6) 

L \rAC J m + 1 \rAC J J 

where — x\g/rAB- Then, given that 

^2(AS) ^ j^j^ p*AP*Bd'xAd'xB^'A^'BSi^AB) , (D7) 

it can be shown that, for two source points in body A and a field point outside the body, 

^2(AA) = - n + —{^1- ^AS^n + Ois) , (D8) 

^ Ua Va 

where y\ = {x — xaY /\^ — Xy^j, and y/i = |x — x/i|, and x\ now denotes the center of baryonic mass of body A. For 
a spherically symmetric body, this agrees with Eqs. (C6) and (C7) of Paper II [T^ . 

Similarly, for a source point and a field point in body A, and the other source point in body B, 

^2(AB) - ^ ^'^ J + 2^2 " 

, ^ik-> f^^^ Hk ^^..tk\ , IruAmB 

-n 



4 r 



3 



<J'=> ( 2XaS'^ - -Yf] + l "^-4ma _fc _ _ ^^^^^ _ ^^J^^k■. 

\ 3 J 4: r-^ 
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+0(s^), (D9) 

where = {x\—Xg)/r, r — jx^— xs|, andx''' — x''—x\; Xa, Ya, and Za denote the superpotentials p*'|x— x'j^d'^x' 
generated by body A only, where p = 1, 3, 5, respectively, and Ia is the moment of inertia tensor of body A. Angular 
brackets <> around indices denote the symmetric trace- free product. We have kept terms up to order in Eq. (ID9P 
because they will ultimately be multiplied by terms that scale as negative powers of s. 

For a source point in body A, a source point in body B and the field point between the two bodies, 

^ij _ TTiAmB ( (ijA ~ ny[yB + ny ^ Sij 



^^(-) = arI) [ AMg) + ) + ^('^^ • ^""''^ 

Other useful identities involving integrals over a sphere surrounding one of the bodies, say body B, include, 

1 d^x 1 1 

-XA • xs - -rAs + 7?.B , (Dlla) 



47r Jo |x-xa||x-xb| 3TZb 2 

i I |x-XA|ix"-XB| = - isk f"^"^ + "^"^^ " ^^"^^ "'^^"^ ■ " i^"^ + "^^^"^^ 

+ (Dllb) 

To illustrate the method used in evaluating J^p*H'^(Px, we consider one integral in Eq. (|D4p . namely 
(47r)^^ Ui^Ui P^^^2)'^^^- Considering first the integral over a sphere of radius TZi surrounding body 1, we expand 
^2(22) powers of x™ — x™ — about the center of mass of body 1. The only term that gives a non-zero result that 
scales as s*^ or lower is the term linear in x. We must therefore evaluate the integral J [/{*^ a;™ d'^x over a sphere 
surrounding body 1. This can be done using Eqs. (jPlip . Combining the result with P2(22) m evaluated at xi using 
Eq. (jD8|l . and keeping terms scaling as or TZ^^ , we obtain the term 

^(nm' - nl'n^ + 2r!f n*^'^) + , ( + Un2n' ~ 2nl'n^ . (D12) 

Ar-^ 307?.ir^ \ r J 

Integrating over a sphere surrounding body 2, we expand the product U^'^U'i about X2 in powers of = x™ — x™. 
But because, inside body 2, ^2(22) scales as s~^, the integral J2 ^2(22)'^'^ ^ scales as s"^~^^, yielding no s^^ or s° 
terms of interest. Finally, to integrate over the domain A4 — Mi — M.2, we note that, because the integral over the 
domain Ai2 yields only positive power scaling, then for our purposes, we can integrate over this domain equally well 
by evaluating the integral over the exterior of the sphere surrounding body 1, J.^ J7{*"' J7{'^P2(22)^'^^' with the external 
potential Ui — rni/r, and the exterior form (jPSp for P2(aA)^ using body 2 as the source. The fact that the chosen 
form of P2^22) singular at X2 does not affect the parts of the integral we are interested in, and the integral is finite. 
The result for this term is 

Combining expressions (|D12p and (|D13[) . we see that the terms proportional to TZi"^ cancel, as they must, leaving the 
result for this term, 

i!!i(r!2«' - ^in^ + 2Ui^n}^^) + i^ff^m* - ^f^ nP + 2r!f . (D14) 

This happens to be antisymmetric on interchange of 1 with 2 (whereby —n^), which is to be expected, since this 
particular term can also be written in the form 

1 f i h i' ui (Pxd^x' ix — x'Y ,„ , 

M 1^ ~ ^ I 



2(47r)2 
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which is manifestly antisymmetric under 1^2. The remaining contributions to Termi in Eq. (|D4[) and Term2 in Eq. 
(|D5p can be evaluated in the same manner. The final result for terms scaling as s^^, together with the point-mass 
contributions, is given by Eq. (|4.6p . 
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